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Excitonic fine structure splitting in quantum dots is closely related to the lateral shape of the 
wave functions. We have studied theoretically the fine structure splitting in InAs quantum dots 
with a type-II confinement imposed by a GaAsSb capping layer. We show that very small values 
of the fine structure splitting comparable with the natural linewidth of the excitonic transitions 
are achievable for realistic quantum dot morphologies despite the structural elongation and the 
piezoelectric field. For example, varying the capping layer thickness allows for a fine tuning of the 
splitting energy. The effect is explained by a strong sensitivity of the hole wave function to the 
morphology of the structure and a mutual compensation of the electron and hole anisotropies. The 
oscillator strength of the excitonic transitions in the studied quantum dots is comparable to those 
with a type-I confinement which makes the dots attractive for quantum communication technology 
as emitters of polarization-entangled photon pairs. 

PACS numbers: 71.35.-y,73.21.La,81.05.Ea 


I. INTRODUCTION 

Excitonic fine structure splitting (ESS) refers to a tiny 
energy splitting of two bright exciton states confined in 
quantum dot (QD) heterostructures with a typical mag¬ 
nitude ranging from units to hundreds /leV. It is mani¬ 
fested in a doublet structure of the exciton recombination 
band. It was observed for the first time in GaAs/AlGaAs 
quantum wells with fluctuating thickness^ and then in 
various QD systems.—^— Soon after its discovery it was 
attributed to the electron-hole exchange interaction^ and 
its finite value was related to the reduced symmetry, 
which needs to be lower than D 2 d^ 

The interest in ESS is triggered by both fundamen¬ 
tal and application point of view. ESS helps to distin¬ 
guish the spectral features originating in the recombina¬ 
tion of exciton (doublet), biexciton (doublet with oppo¬ 
site polarization-energy dependence), and trion (singlet). 
Benson’s proposal of the source of entangled photon pairs 
relying on zero ESSi has called for the preparation of 
QD systems with low ESS. Using (111) substrates for 
the growth of InAs QDs reduced both structural asym¬ 
metry and piezoelectric contributioni^ Another attempt 
involved strain-free GaAs/AlGaAs QDs with zero piezo¬ 
electric field, which, however, still exhibited a finite ESS 
due to structural elongationi^ Post-growth annealingiS 
of InAs QDs allowed to decrease ESS from 96 pieV to 
mere 6 /xeV. Another class of approaches is based on 
in-operation tuning, where the originally large value of 
ESS is reduced by applying the external field: elec- 


tricfiid^ magnetic^d^ or strain. The external strain 
field allowed to reach ESS below experimental resolution 
in GaAs/AlGaAs QDs j^^d^ the simultaneous application 
of electric field allowed for a more powerful symmetry 
restoration and rather universal recovery of low ESS.— 

Various effects contributing to the ESS can be di¬ 
vided into two classes based on the involved length scale: 
atomic and macroscopic. Atomic-scale effects are con¬ 
nected with the irregularities of the crystal lattice such 
as the interfaces, particular elements distribution in al- 
loySfi^ charged defects^^ etc. The magnitude of these 
effects is still subject of investigation; the atomistic simu¬ 
lations based on the tight-binding methodic predict con¬ 
siderably larger values than those relying on the empir¬ 
ical pseudopotential method— In general, atomic-scale 
effects are weak compared to those on macroscale. For 
example, the magnitudes of about 1 fj,eV are reported 
for a specific alloy distribution in the AlGaAs barrier— 
of GaAs QDs. The effect is more pronounced when 
the dot material is an alloy, which should be therefore 
avoided when aiming at low ESS. A lower bound of sev¬ 
eral /xeV was predicted for strain-tuned ESS in ternary 
Ino.6Gao.4As QDs— By macroscopic scale we mean for 
the purpose of the foregoing discussion that the charac¬ 
teristic length of the effect is comparable with the dimen¬ 
sions of a QD and the underlying crystal lattice is per¬ 
ceived as a homogeneous environment. Thus, the crystal 
symmetry is no longer relevant and the finite values of 
ESS are now related to the symmetry lower than C 4 , 
i.e., to the lateral elongation of the wave functions (e.g. 
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envelope functions of the k • p theory). Principal contri¬ 
butions to the FSS on macroscopic scale arise from the 
asymmetric (elongated) shape of a QD and piezoelectric 
field. Further, it is possible to use the external strain 
field to induce the anisotropic effective mass tensor and 
modify the elongation of the hole wave functions and the 
related value of FSSf^^ 

Further information is contained in the polarization 
properties of the exciton doublet. Simple considerations 
assuming a purely heavy-hole exciton in an elliptic-disk¬ 
shaped QD^ predicted a linear polarization of both tran¬ 
sitions with the low-energy component polarized paral¬ 
lel with the long QD axis and the high-energy compo¬ 
nent having the orthogonal polarization. Typically, both 
structural-elongation and piezoelectric axes are parallel 
with the crystal axes [110] and [lIO], and so are the po¬ 
larizations of both components. However, in some struc¬ 
tures with shallow irregular confinement potential, such 
as quantum well thickness fluctuations, stochastic po¬ 
larization directions were observedi^ Further, when the 
light-hole contribution to the exciton ground state be¬ 
comes important, the polarization orthogonality of both 
components is lost i^^’^^ 

We focus here on QDs with type-II conhnement, in 
which one type of charge carriers is confined in QD vol¬ 
ume and the other in the barrier close to the QD vicin¬ 
ity. The particular system of interest are InAs QDs with 
a thin GaAsi_ySby overlayer embedded in GaAs. One 
reason for selecting this material system is the possibility 
to induce a smooth crossover between type-I and type-11 
conhnement regime simply by changing y, the crossover 
values between 0.14 and 0.18 have been reportedi^^— 
The other is that it belongs to the minority of systems 
with holes bound outside. Owing to their large effective 
mass the holes are more susceptible to the local potential 
prohle or external perturbations, offering a larger poten¬ 
tial for tuning their wave functions and the related FSS. 
The photoluminescence of GaAsi_ySby capped QDs is 
rather intense despite the type-II conhnement with the 
radiative lifetimes as low as 10 nsi^ The strain-reducing 
effect of GaAsi-ySby layer together with the surfacting 
effect of antimony allow to increase the emission wave¬ 
lengths of standard InAs QDs and reach the telecom¬ 
munication wavelength of 1.3 and 1.55 Various 

shapes of GaAsi_ySby QDs have been reported, includ¬ 
ing a lens^E or a pyramid with a graded In concentra- 
tioni^S Notably, the hole wave function is expected to be 
composed of two segments localized in the minima of the 
piezoelectric potential<^ 

In this work we present a theoretical study of excitonic 
hne structure splitting of InAs QDs with GaAsi_ySby 
overlayer. We propose a method to tune the FSS by set¬ 
ting the thickness of the GaAsi_ySby layer. The values 
comparable with the natural linewidth can be achieved 
even in low-symmetry QDs. The paper is organized as 
follows: In Section El a theory of FSS is described. To 
gain a qualitative understanding of the relations between 
the wave functions and FSS we discuss in Section unia 


simplified single band model with Gaussian wave func¬ 
tions. The full calculations are presented in Section HVl 
We summarize and conclude in Section IVl 


II. THEORY 

The single particle states were calculated within the 
eight-band k • p theoryj^LS^ in which the wave functions 
are expanded into products of periodic parts of Bloch 
functions rtf, in the F point and corresponding envelope 
functions Xb, 

V'(i’) = W6WXfe(r). (1) 

In this equation b is the band index, the bands x, y, z 
correspond to the valence band Bloch waves which are 
antisymmetric with respect to the corresponding mirror 
plane and s corresponds to the conduction band Bloch 
wave. Following usual conventions, z denotes the growth 
direction. The calculations include the effects of the elas¬ 
tic strain via the Pikus-Bir Hamiltonian^ and the piezo¬ 
electric field. The numeric simulations were performed 
with Nextnano 3D^ which employs the finite difference 
method. The simulation space was discretized with a 
step of 1 nm. 

Once the single particle states are calculated, it is con¬ 
venient to use them as a basis for the exciton state |A1). 
First the Slater determinants \X{ci,vj)) = clyyjlO) are 
formed, where |0) is Fermi vacuum state (empty quan¬ 
tum dot), creates an electron in ith conduction state 
and Cyj annihilates an electron in jth valence state, the 
corresponding single-particle wave functions are denoted 
'4’ci, '>Pvj, respectively. For the calculations of FSS we 
used four Slater determinants formed from the ground 
hole and electron states. Following Ref. i (Eq. 2.3) the 
Hamiltonian matrix elements read 

{X{ci, vj)\H\X{ck, vl)) = 

{Ei — Ej)SikSji + C{ci, vj, ck, vl) + EX{ci, vj, ck, vl) 

where Ei is the energy of ith single particle state, C 
represents the direct Coulomb interaction, and EX rep¬ 
resents the exchange interaction. 

Defining 

5'c.«(r) = Xc,b(j')x«,b(r) 

be{s,a:,y,z}®{t,4.} 

and vector T with the components 

Tc,y(r) = ^ [x*c,sa{^)Xv,xa{r)+X*c,xa{^)Xv,s^{r)] 

and T^, defined analogously {Eg is the fundamental 
band gap and P is one of the Kane’s parameter related to 


the non-vanishing coordinate matrix elements {a;|x|s) = 
P/Eg), we can write 


C{cl, vl] c2, v2) 



-—X 
47re 


X 


dri 


dr2-. -r5'ci,c2(i’i)5't,2,t,i(r2) 

|ri - r 2 | 


(e denotes the elementary charge and e the dielectric 
function). The exchange Coulomb interaction term EX 
can be expressed as a sum of the following three terms: 


EXo{cl, ul; c2, v2) = 




— X 
47r£r 


dr2-. -rS'ci,„i(ri)5'„2,c2(i’2), 

|ri2| 


( 2 ) 


EX^{cl,vl-c2,v2) = ^ I dr, j dr 2 j^x 
X ri2 • [<S'.u2,c2(r2)Tcl,^l(ri) — Scl,vl{'l‘^l)T^v2,c2{^2)] 


( 3 ) 


and 


i?X2(cl,'!;l;c2,?;2) =/ dri / dr2 , ^ x 

dTrery J |ri2p 

X '^cilii^i)Tv 2 l 2 i^ 2 )[Sc,fi\ri 2 \‘^ - (4) 

a,^^x,y,z 

where ri 2 = ri — r 2 and Sap is the Kronecker delta. We 
note that S is non-zero only when the mixing of valence 
and conduction bands is taken into account. Thus, only 
the third term of the multipole expansion, Eq. |4l con¬ 
tributes to the FSS when this mixing is neglected, e.g. 
when single band or six-band^ k.p theory is used to ob¬ 
tain the wave function of individual particles. However, 
as the scaling of the terms with the linear extension of 
the wave function L goes as EXq ~ 1/L, EX, ^ 1/T^, 
EX 2 ^ 1/T^, the low-order terms are important in par¬ 
ticular in larger QDs. 


III. MODEL OF GAUSSIAN WAVE FUNCTIONS 


Before treating realistic quantum dots with the full- 
complexity model, it is worth to provide an intuitive 
understanding of the relation between the topology of 
the excitonic wave function and the value of FSS. To 
this end we employed a simplified model with the exci- 
ton composed of a single Slater determinant, neglected 
band mixing, and the electron and hole densities having 
the form of three-dimensional Gaussian functions. The 
electron and hole envelope functions read 


Xe,h{r) oc exp 


(x-xo)^ (y-yo)'^ (z - zo) 
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FIG. 1: (color online) FSS as a function of the lateral 
elongation of the wave functions, dehned as E™ = L^/Ly. 
(orange line) A smaller dot with the extension parameters 
Lx X Ly = 25 nm^, Lz = 2 nm. The lateral extensions are 
varied preserving the value of their product (and thus the 
volume and vertical aspect ratio of a QD). Electron and hole 
wave functions are elongated equally, (green line) A larger 
dot with the extensions twice larger than for the smaller dot, 
i.e., LxXLy = 100 nm^, Lz = S nm. (magenta line) A smaller 
dot. Lx X Ly = 25 nm^, Lz = 2 nm, with opposite electron 
and hole elongation. For electrons, Eu, has a constant value 
of 2/3. For holes, E^j is varied and FSS is displayed as a func¬ 
tion of the hole elongation, (insets) The insets schematically 
depict the topology of the wave functions. Hatched magenta 
ellipses correspond to the electron wave functions. 


where Lx,y,z determine the spatial extensions of the wave 
functions and Xo,yo, zq are the coordinates of the particle 
barycenter. As the band mixing is neglected, only the 
dipole-dipole exchange term (Eq. |3|) contributes to the 
total FSS. 

A crucial parameter for FSS is the lateral elongation of 
the envelope functions defined as Ey, = Lx/Ly. It follows 
directly from Eq.[5]that for non-elongated envelope func¬ 
tions {Lj: = Ly) FSS acquires a zero value. The depen¬ 
dence of the FSS on the elongation is shown in Figure [T] 
In order to isolate the effect of the elongation and avoid 
unintentional variation of other parameters, we preserved 
the volume and the effective vertical aspect ratio of the 
model dots, i.e, values of L^ and the product L^ x Ly were 
kept constant. First, we assumed the same envelope func¬ 
tion for both electron and holes (orange and green line). 
Such case corresponds e.g. to strain-free GaAs/AlGaAs 
dotsi^ FSS exhibits a monotonically increasing concave 
dependence on the lateral elongation Ey,. To demon¬ 
strate the effect of the QD volume, we show FSS for a 
smaller dot (extension parameters Ly, x Ly = 25 nm^, 
Lz = 2 nm) and a larger dot with two-times larger di¬ 
mensions (e.g., eight-times larger volume). The values 
of FSS for a larger QD are exactly eight-times smaller. 
The inverse proportionality of FSS to the QD volume or 
to the third power of a characteristic linear dimension L 















can be directly inferred from Eq. |4l With the band mix¬ 
ing taken into account, additional terms proportional to 
l/L and 1/L^ emerge. However, the 1/L^ or \/V scal¬ 
ing law (V representing a volume of the QD) has been 
recently demonstrated experimentally in realistic strain- 
free GaAs/AlGaAs QDs.— FSS values exceed the natural 
linewidth of the exciton recombination lines (up to units 
of fJ.eV) even for a modest elongation. For example, for 
= 1.2 we predict FSS of 11 /reV (1.4 /xeV) in the 
smaller (larger) QD. We note that for the QDs studied 
in Ref. we found the values of x Ly between 11 
and 36 nm^ and between 1 and 2 nm. The smaller 
dot case thus corresponds well to realistic GaAs/AlGaAS 
QDs. 

Next, we introduce an important concept of the com¬ 
pensated elongation. For a suitable exciton wave func¬ 
tion topology, FSS can attain the zero value even in the 
system that lacks the required symmetry C 4 v. We con¬ 
sider the electron envelope function to be elongated in the 
direction perpendicular to the elongation of the hole en¬ 
velope function, E^je = 2/3 (the subscripts e, h are used, 
when required, to distinguish the parameters of electrons 
and holes, respectively). The extension parameters corre¬ 
spond to the smaller dot: x Ly = 25 nm^, = 2 nm. 

FSS is plot as a function of the hole elongation E^jh > 1 
in Fig. [T] (magenta line). The prominent feature of the de¬ 
pendence is the zero-value minimum at E^h = 3/2, (i.e., 
the inverse of the hole elongation. Intuitively, this can be 
described as the mutual compensation of both electron 
and hole elongations. The integral in Eq.|4] attains a zero 
value, which is however not related to the symmetry. In 
realistic QDs, the condition of the inverse elongation does 
not hold (due to band mixing or different volume of the 
envelope function of electrons and holes) but the effect 
is preserved. The minimum value of FSS can be larger 
than zero in case of QDs with the irregular shape. The 
effect of the compensated elongation has been already 
demonstrated experimentally utilizing the anisotropic ex¬ 
ternal strain to vary the elongation of the hole envelope 
functionJ^ In the following we will demonstrate that the 
elongation of the hole envelope function can be efficiently 
varied in type-II InAs QDs with GaAsSb capping layer. 

The transition between type-I and type-II confinement 
in GaAsi-ySby capped QDs is accompanied by the split¬ 
ting of the hole wave function into two segments, evenly 
spread along the central electron wave function in QDs 
with a sufficient symmetry (Fig. [6]3,c). The behavior of 
FSS under such transition is shown in Fig. [2] When the 
connecting line of the segments is parallel to the elonga¬ 
tion axis (orange line), small segment shifts from the cen¬ 
tral position effectively enhance the hole elongation and 
consequently the FSS. For larger shifts the wave func¬ 
tion disintegrates; now the effect of increased distance 
of electron and hole prevails resulting into a decrease of 
FSS. The same behavior is predicted for non-elongated 
wave functions (green line), where the FSS dependence 
starts at a zero value, increases as the holes become ef¬ 
fectively elongated, and decreases when the separation 



FIG. 2: (color online) FSS for a hole wave function composed 
of two segments. The extensions of both electron and hole 
wave functions fulfill L^x Ly = 25 nm^, Lz = 2 nm, with the 
lateral elongation of 5/4 (orange line, squares), 1 (green 
line, circles), 5/6 (magenta line, diamonds), and 5/8 (maroon 
line, triangles). The electron wave function is composed of 
a single Gaussian while the hole wave function is gradually 
split along x into two identical segments evenly positioned 
around the central electron wave function. FSS is plotted 
as a function of the distance between the barycenters of the 
segments. Insets schematically depict the position and shape 
of the segments. 


effects prevail. When the connecting line of the segments 
and the elongation axis are perpendicular, the effective 
elongation of the holes decreases as the segments are sep¬ 
arated, and so happens with FSS. Depending on the mag¬ 
nitude of the original elongation, two possibilities exist: 
(1) The hole eventually becomes elongated in opposite di¬ 
rection (magenta line; note the magenta insets of Fig. [2] 
schematically depicting the change in the elongation di¬ 
rection). FSS goes through a zero value and starts to 
increase again. Finally, the separation effects prevail and 
FSS decreases. (2) When the original elongation is large, 
the hole disintegrates before the elongation direction is 
changed (maroon line). In such case a monotonously de¬ 
creasing dependence of FSS on the segment distance is 
observed, governed first by the decrease of the effective 
elongation and then by the separation effects. 

In strongly asymmetric QDs the minima in confine¬ 
ment potential corresponding to both segments can dif¬ 
fer considerably and a single-segment hole wave function 
displaced from the central electron wave function can 
be formed. In such case there is no effective change 
in the elongation and the separation effect leads to a 
monotonous decrease of FSS as in the case of large per¬ 
pendicular elongation (not shown). 
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IV. REALISTIC QUANTUM DOTS 

We will now focus on realistic InAs QDs with a 
GaAsi-ySby capping layer. We will show that FSS in 
such structures can be tuned by the thickness of the 
GaAsi-ySby layer, allowing its decrease below the natu¬ 
ral linewidth of the exciton transitions. 

To show the universality of the tuning approach, 
three QD geometries will be considered here, denoted as 
pyramidal, symmetric lens-shaped, and elongated lens¬ 
shaped. The pyramidal QD is adopted from Ref. [s^ 
and has the shape of a pyramid with the base length of 
22 nm, height of 8 nm, and the trumpet indium composi¬ 
tion profile within the pyramid. For the other structures 
we assume QDs composed of pure InAs. The symmet¬ 
ric lens-shaped QD is modeled as a top of a sphere with 
the base radius of 8 nm and the height of 4 nm. The 
prominent cause of the lateral asymmetry and contribu¬ 
tor to FSS of many QD systems is a structural elonga¬ 
tion. So far no elongation was reported for InAs QDs 
with GaAsi-ySby overlayer, which is in striking contrast 
with InAs QDs capped by pure GaAsj ^^i^^ This can be 
attributed to the surfacting effect of antimony but it is 
also possible that the elongation has been overlooked as 
the methods involved in experimental studies were insen¬ 
sitive to it. Therefore, we consider in our study also the 
possibility that QDs are elongated. In accordance with 
GaAs capped InAs QDs we select the direction [110] as 
the main elongation axis and quantify the elongation by 
the ratio of characteristic lateral dimensions along [110] 
and [110], denoted as Eg in the following (the subscript s 
is used to differentiate the structural elongation from the 
wave-function elongation used in the previous section). 
The elongated lens-shaped QDs are formed from the lens¬ 
shaped dot by its stretching along [llO] and compressing 
along [110] by the same factor so that the QD volume 
and height are preserved. All QDs are capped with the 
GaAso, 8 Sbo .2 layer of a certain thickness and further em¬ 
bedded in GaAs. 

The topology of the wave functions is closely connected 
with the effective confinement potential, which is con¬ 
tributed by the band-edge offsets, strain field, and piezo¬ 
electric potential. We construct the potential from the 
eigenvalues of the pointwise diagonalized Hamiltonian 
(terms containing the spatial derivatives are discarded 
and the Hamiltonian is then diagonalized at each point 
of the simulation grid) so the strain-induced band mix¬ 
ing is already involved. The hole potentials in all type-H 
QDs discussed further in the paper exhibit qualitatively 
similar features. We will present them for the exemplar 
lens-shaped QD with the GaAsSb layer thickness of 5 
nm and with the elongation Eg of either 1 or 2. Figure [3] 
shows the potential profile without the piezoelectric con¬ 
tribution. The potential is given from the electron view, 
the largest values correspond to the minima of the hole 
confinement. Zero energy is set to the valence band edge 
of bulk unstrained InAs. Two local minima of the hole 
confinement potential are formed in the GaAsSb layer 




FIG. 3: (color online) Effective confinement potential for 
holes (colorbar scale in meV) without the contribution of the 
piezoelectric field for the lens-shaped QD. (a) Eg = 1, plane 
(lIO) through the QD center, (b) Eg = 2, plane (001) just 
above the QD base. Boundaries between different materials 
are schematically depicted. The potential is given from the 
electron view; the holes are confined near the largest values 
represented by the red/orange spots. 



FIG. 4: (color online) Piezoelectric potential (colorbar scale 
in meV) for the symmetric lens-shaped QD. (a) plane (110) 
through the QD center, (b) plane (001) just above the top of 
the QD (i.e., above the horizontal nodal plane). 


along the sides of the QD and above its top [Fig. [SKa)]. 
The top-minimum is of the light-hole character and there¬ 
fore penalized by the quantum confinement. The ground 
hole state will be localized in the side-minimum which 
forms a ring around the QD [Fig. EKb)], in which weak 
variations of the potential are present with two rather 
shallow absolute minima along the long QD side (i.e., in 
[110] direction from the QD center). Depending on the 
magnitude of these variations, the wave function might 
form a ring or be split into two segments. 

The piezoelectric field has an octopole shape shown 
in Figure S] Its contribution is rather important as its 
magnitude of about 50 meV is comparable to the varia¬ 
tions of the rest of the confinement potential inside the 
GaAsi-ySbylayer. The horizontal nodal plane of the 
piezoelectric octopole lies close to the side-minimum of 
the confinement potential. The piezoelectric potential 
therefore tends to split the wave function of the holes in 
the side-minimum into two segments situated along [110] 
below the nodal plane or along [110] above the nodal 
plane. 

The total confinement potential is shown in Figure [Sj 
The following morphologies of the hole wave function are 
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FIG. 5: (color online) Total confinement potential for holes 
(colorbar scale in meV) including the contribution of the 
piezoelectric field for the lens-shaped QD. (a,c) Plane (001) 
below the nodal plane of the piezoelectric potential (just 
above the QD base), Eg = 1 (a) and 2 (c). (b,d) Plane (001) 
above the nodal plane of the piezoelectric potential (3 nm 
above the QD base), Eg = 1 (b) and 2 (d). Black circles and 
ellipses display the QD boundary at respective heights. The 
potential is given from the electron view; the holes are con¬ 
fined near the largest values represented by the red/orange 
spots. 


FIG. 6: (color online) (a,b,c) Planar probability density of 
the hole ground state in the lens-shaped QD (height 4 nm, 
radius 8 nm, symmetric), plane (001) (integrated over [001]) 
for the GaAso. 8 Sbo ,2 layer thickness of (a) 2 nm, (b) 4 nm, and 
(c) 8 cm. The QD boundary at the base-height is displayed 
with the white line, (d) FSS (thick solid lines with symbols, 
left) and transition probability (thin dashed lines) between 
the ground electron and hole state in the symmetric lens¬ 
shaped QD (orange) and the elongated lens-shaped QD with 
Eg = 1.1 (green) as functions of the GaAsSb layer thickness. 
The red labels denote the points corresponding to the maps 
of the probability density. 


possible: (i) inside a QD (type I), (ii) a ring-like shape 
along the QD when the piezoelectric field is too weak 
to localize the holes in its minima, (iii) two segments at 
the dot base situated along [110] [Fig. [5l[a,c)], (iv) two 
segments at the QD sides above the piezoelectric nodal 
plane situated along [110] due to the piezoelectric field 
[Fig-EIb)] or along [110] when the structural elongation 
along [lIO] prevails [Fig. [5l[d) is close to that case]. In 
pyramidal QDs with the trumpet In composition pro¬ 
file, both the side-minimum of the confinement potential 
and the piezoelectric octopole are shifted up towards the 
region of large In content. In short, there is a rich vari¬ 
ety of the hole wave function morphologies. Variation of 
the parameters such as the thickness or the composition 
of the GaAsSb layer are supposed to induce transitions 
between those morphologies. For example, switching be¬ 
tween the deep narrow minima below the nodal plane 
[Fig. [5][a,c)] and shallow broad minima above the nodal 
plane [Fig. [5l[b,d)] shall be achievable. Considering the 
effect of the compensated elongation, this opens an in¬ 
teresting prospect for the tuning of FSS. 

Figure |6][d) shows the dependence of FSS on the thick¬ 
ness of the GaAsi-ySby layer in two lens-shaped QDs: 
symmetric and weakly elongated {Eg = 1.1). We will 
first discuss the case of the symmetric QD. The electron 


wave function is weakly elongated in [110] direction and 
as it resides within the QD, its variation with the thick¬ 
ness of the GaAsi-ySby layer is negligible. For a thin 
GaAsi-ySby layer (up to 3 nm) the ground hole wave 
function resides inside the QD, too. It experiences the 
bottom part of the piezoelectric octopole and is thus elon¬ 
gated in [110], as shown in Fig. |6l[a). The polarization 
of the lower exciton component is [110]. With increasing 
thickness a hole ground state gradually shifts into the 
GaAsi-ySby layer and also slightly upwards (for about 
1.3 nm for the full range of thicknesses). Consequently, 
it becomes split by the upper part of the piezoelectric oc¬ 
topole into two segments along [110]. For the thickness 
interval between 3 and 5 nm, the segmented wave func¬ 
tion behaves as effectively elongated in [lIO] [Fig. [SKb)] 
and the lower exciton component is also polarized along 
[110]. For the thickness values above 5 nm, the segments 
are well separated and the elongation of each segment 
in [110] determines the overall symmetry of the wave 
function. The lower exciton component is again polar¬ 
ized along [110]. Thus, we distinguish three regions of 
different exciton polarization. At each of the two transi¬ 
tions between those regions, exciton levels cross and FSS 
is reduced to zero. We note that minimum values ob¬ 
tained in our calculations are non-zero due to the hnite 
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FIG. 7: (color online) FSS (thick solid line with symbols, 
left) and transition probability (thin dashed line) between the 
ground electron and hole state in the pyramid-shaped QD as 
functions of the GaAsSb layer thickness. 

step in the thickness dependence and read 0.1 ^eV and 
1 ^eV for the first and second transition, respectively. 
The observed behavior of FSS corresponds well to the 
qualitative prediction of the model of Gaussian function 
for LxjLy = 5/6 (cmp. Fig. O magenta line). 

Similar behavior is observed in a lens-shaped QD 
weakly elongated in [110] [Figure [SKd), green lines]. The 
first region is now missing, as the elongation of the 
hole wave function within the QD volume (for a thin 
GaAsi-ySby layer) is now [110] because the structural 
elongation dominates over the piezoelectric contribution. 
The zero FSS is reached for a GaAsi_ySby layer thick¬ 
ness of 7 nm. Figure |7| shows the FSS dependence on the 
thickness of the SRL layer in a pyramidal QD with the 
trumpet profile (height of 8 nm and base size of 22 nm). 
As in the case of the lens-shaped QDs, the hole function 
changes its effective elongation during the crossover be¬ 
tween the type-I and type II confinement and in turn, ex- 
citon levels cross and FSS is reduced towards zero. Thus, 
variations of the GaAsi_ySby layer thickness present a 
universal approach to tune and reduce FSS. Further, its 
tuning power is considerable as it can compensate even a 
structural elongation, as was demonstrated for the elon¬ 
gated lens-shaped QD. For a typical lifetime of InAs 
QDs of 1 ns, the natural linewidth of the exciton re¬ 
combination is 4 /reV; experimentally predicted spectral 
linewidth approaches 100 /xeV4£ Both values are larger 
than the minimum accessible values predicted here for 
the GaAsi-ySby capped QDs. 

Finally we argue that GaAsi_ySby capped InAs QDs 
are good emitters even in type-II confinement regime, 
contrasting the other type-II material systems in which 
the studies of the optical properties of individual QDs 


are rare and challenging i^^d^ There is an experimen¬ 
tal evidence of intense photoluminescence from the 
GaAsi-ySby capped InAs QDs,— about 4-times weaker 
than for the reference type-I sample and about 3-times 
stronger after rapid thermal annealing. Roughly 15-times 
higher exciton lifetime has been reported (0.7 ns in type 
I compared to 11.2 ns in type II). We support these 
observations with the calculated transition probabilities 
(square moduli of matrix elements) shown in Figs. |6][d) 
and [71 The values are normalized to the full overlap of 
the envelope functions. Typical values for type I QDs 
are in the range 0.5-0.7. The values corresponding to 
the minima of FSS are 0.37 and 0.15 for the symmetric 
lens-shaped QD, 0.10 for the elongated lens-shaped QD, 
and 0.12 for the pyramidal QD. Thus, roughly 5-times 
weaker photoluminescence as compared to type-I QDs 
is expected. For this reason we are convinced that the 
extraction of individual photons or photon pairs from in¬ 
dividual GaAsi-ySby capped QDs will be experimentally 
feasible. 


V. CONCLUSIONS 

In GaAsi-ySby capped InAs QDs, lateral symmetry 
of the hole wave functions can be to a large extent influ¬ 
enced by the thickness of the GaAsi_ySby layer. In par¬ 
ticular, during the crossover between type-I and type-II 
confinement regimes, the hole wave function is shifted up¬ 
wards across the nodal plane of the piezoelectric octopole, 
which is accompanied by the change of the direction of 
lateral elongation. Due to the mechanism of compen¬ 
sated elongation, a crossing of the bright exciton levels 
and a reduction of FSS to zero is predicted for certain 
thicknesses of the GaAsi_ySby layer. Low natural FSS 
and efficient photoluminescence make the GaAsi_ySby 
capped InAs QDs attractive as a possible source of en¬ 
tangled photon pairs. 
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